Hidden Quasiparticles and Incoherent Photoemission Spectra in Na2lr03 
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We study two Heisenberg-Kitaev t-J-like models on a honeycomb lattice, focusing on the zigzag 
magnetic phase of Na2lr03, and investigate hole motion by exact diagonalization and variational 
methods. The spectral functions are quite distinct from those of cuprates and are dominated by 
large incoherent spectral weight at high energy, almost independent of the microscopic parameters 
— a universal and generic feature for zigzag magnetic correlations. We explain why quasiparticles 
at low energy are strongly suppressed in the photoemission spectra and determine an analog of 
a pseudogap. We point out that the qualitative features of the predominantly incoherent spectra 
obtained within the two different models for the zigzag phase are similar, and they have remarkable 
similarity to recently reported angular resolved photoemission spectra for Na2lr03. 

PACS numbers: 75.10.Jm, 72.10.Di, 75.25.Dk, 79.60.-i 
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Frustrated spin systems have long served as a rich 
source of exotic phenomena such as unconventional mag- 
netic order or spin liquid behavior fH, f3| ■ A beautiful re- 
alization of a spin liquid is found in the Kitaev model [3] , 
where bond-dependent Ising interactions lead to strong 
frustration on the two-dimensional (2D) honeycomb lat- 
tice. Here antiferromagnetic (AF) or ferromagnetic (FM) 
couplings are equivalent; a FM realization was originally 
suggested in Ref. f3|. The Kitaev model is exactly solv- 
able and has spin correlations that vanish beyond nearest 
neighbor (NN) spins [5] . This spin liquid is robust against 
weak time-reversal invariant perturbations, including the 
Heisenberg ones [6], similar to the spin-orbital liquid in 
the SU(4) symmetric Kugel-Khomskii model Q, but in 
striking contrast to the 2D compass model 

Such peculiar, highly anisotropic interactions are be- 
lieved to be realized in A2h:0^ {A =Na,Li) iridates where 
strong spin-orbit coupling generates Kramers doublets of 
isospin states |^ [ld| ■ These systems are Mott insulators 
as confirmed by the electronic structure [ll, 12|. It has 
been suggested that effective S = \ pseudospins, stand- 
ing for locally spin-orbital-entangled t2g states l3, 1^ . 
interact via competing AF Heisenberg and FM Kitaev ex- 
change couplings between NNs in the Heisenberg-Kitaev 
(HK) model — this scenario is consistent with the mag- 
netic susceptibility fTBI and with resonant inelastic x-ray 
scattering (RIXS) [16]. Surprisingly, the theoretical pre- 
dictions of a spin liquid or stripy AF phase Q were not 
confirmed but instead zigzag AF order, consisting of stag- 
gered FM zigzag chains, was found in Na2lr03 [17-20j. 
To stabilize this phase next-nearest neighbor (NNN) and 
third nearest neighbor (3NN) AF Heisenberg terms [2l| 
have been invoked in the HK model [22J- Recently a sim- 
pler scenario, including e„ orbitals in the exchange paths, 
has been proposed 23|, |2J| ; there, the zigzag AF phase 



emerges in a broad range of parameters in the HK model 
when the signs of all NN exchange terms are inverted. 

In this Letter we investigate whether hole motion in 
an effective HK model can explain recent angle-resolved 
photoemission spectroscopy (ARPES) experiments for 
Na2lr03 [ll|. One expects that a hole might move co- 
herently along the FM zigzag chain [Fig. 1(a)], simi- 
lar to free hole propagation in states with ferro-order of 
t2g orbitals [25|]. It is therefore quite surprising that the 
ARPES spectra for Na2lr03 are dominated instead by 
incoherent spectral weight found predominantly at high 
energy jTlj. This poses several intriguing questions: (i) 
Can one distinguish by ARPES the two zigzag states re- 
alized by qualitatively different interactions? (ii) Is there 
any similarity between the present ARPES results and 
those known for cuprates [26], described by quasiparti- 
cles (QPs) within the t-J model |27l.[28|. or are the spec- 
tra dominated by entangled spin-orbital excitations (29|? 





Figure 1. (color online). Honeycomb lattice of Na2lr03: (a) 
cluster of A'^ = 24 sites, and the elementary translations ei(2) 
within two sublattices A or B (for two NN A atoms a = 
1). Exchange couplings are: in Heisenberg {Jh) and Kitaev 
{Jk) (solid); Heisenberg NNN (J2) and 3NN (J3) (dashed), 
(b) First Brillouin zone with high symmetry points: F, My = 
(tt, — 7r/\/3), Kz ~ (47r/3, 0); in exact diagonalization M-y and 
are equivalent (7 G {x, y, z}) and called M and K. 
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(iii) Is the incoherent hole scattering on spin excitations 
as important here as between the FM planes of LaMnOs 
(sot? We answer them below and show that weak QPs 
exist but are hidden for the present honeycomb lattice so 
that ARPES reveals the incoherent processes. 
We consider the following t-J-like model (< > 0), 



+ JhJ2s,-S, + JkJ2S^S]. (1) 



on the honeycomb lattice [Fig. [IJa)], called model I. The 
hopping "Ht of composite fermions with pseudospin flavor 
a \3l\ (contributions from direct d-d and d-p-d hopping 
[32| ) occurs in the restricted space without double occu- 
pancies due to large on-site Coulomb repulsion U, and 
the exchange Hj stands for the HK model — with FM 
Heisenberg {Jh < 0) and AF Kitaev {Jk > 0) exchange, 



Jh 



-J(l-a), JK = 2Ja. 



(2) 



Here J is the energy unit and a E [0, 1] interpolates be- 
tween the Heisenberg and Kitaev exchange couplings for 
NN i spins {Si}; the zigzag AF phase (Fig. 2) is found in 
a broad range of a jlJ] • We also investigate the spectral 
functions in model II, where exchange couplings — are 
extended by NNN and 3NN terms J„ = J(l - a)j„ > 0: 



{ij}eNNN {ij}e3NN 



SrS,]. (3) 



Although our aim is to present the spectral functions 
obtained by exact diagonalization (ED) for a hole in 
the quantum-fluctuating zigzag AF phase, we begin with 
describing the physical processes and resulting spectral 
properties of a hole inserted into a fully-polarized ground 
state |0), see Fig. 2(a). The hole hopping is isotropic, 
but for convenience we distinguish intrachain hopping t 
and interchain hopping t±. Free hole propagation, see 
Fig. 2(b), occurs along the ID FM zigzag chain when 
t± — 0. It involves two types of Ir sites which belong to 
sublattices A and B, see Fig. 1(a). We are interested in 
the spectral properties measured in the ARPES exper- 
iment with the hole creation operator ct , and we also 

consider the hole creation on sublattice A, dt 1331 1: 

fct — 



As we shall see, the sublattice aspect is rather subtle 
and responsible for the hidden QP states in the ARPES 
experiments for the zigzag phase. The spectral functions 
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Figure 2. (color online). Artist's view of hole propagation in 
the zigzag phase of Na2lr03 (arrows), (a) The hole (circle) 
may either (b) propagate along the FM chain by intrachain 
hopping t, or (c) create a magnon by interchain hopping t±. 
Either this spin defect (d) or the hole (e) can move along its 
chain, (f) After two steps (d) + (e) the hole and the defect 
may recombine. Broken bonds are indicated by dashed lines. 



correspond to the physical Green's function Gc{k,uj) for 
/ = c, or to the sublattice Green's function Gd{k,Ld) for 
f = d here Eq is the energy of the ground state |0). 

Below we describe and use a variational approach well 
adapted to the perturbative regime {\t\ <C J) to gain 
more physical insights. We split the HK t-J model Eq. 
(1) into an exactly solvable part "Ho and the perturbation 
V; Htj ^ Ho + V. When exchange interactions are Ising- 
like, either in spin [s?] or in orbital [25| systems, the num- 
ber of magnons (orbitons) is conserved and the classical 
ground state |0) is exact. Here we use the same strat- 
egy to develop a variational treatment [sB'l and include 
in Hq all terms that conserve the number of magnons, 
while V includes the interchain hole hopping which cre- 
ates or removes a spin excitation [Fig. 2(c)] together with 
exchange terms generating magnon pairs. 

When a hole moves to a neighboring chain, it creates 
a spin defect with energy = \Jh\, see Fig. 2(c). The 
defect propagates as a magnon [Fig. 2(d)], and the hole 
can also move by hopping t [Fig. 2(e)] along the FM 
chain. Both processes generate two parallel spins on a 
vertical (interchain) bond and increase the magnon en- 
ergy to Em = \Jh\ + -^{Jk + Jh), but a hole and a spin de- 
fect may also meet at one vertical bond [Fig. 2(f)], which 
decreases the magnon energy back to £m'' ■ All these pro- 
cesses cause incoherence. Other states with several spin 
defects cost too much energy when J ^ t — these states 
are neglected in the one-magnon variational treatment. 

In the zigzag phase we write the Dyson's equation, 
G(w) = G°{uj) + G(w)VG°(a;), for the resolvent, G{uj) = 
{ut — irj + Eo — T-L}^^. The Green's functions Gd(k, w) ob- 
tained from it are 2x2 matrices for two sublattices, A and 
B. The unperturbed Green's function, G°(k, w), is found 
exactly — it describes free propagation of a hole along 
the FM chains, and depends on ID momentum kx [Fig. 
1(b)], e^_|_ — eoT 2tcos^kx, with eq = jJ{l + a) being 
the energy of a static hole. The two ± states for each kx 



Figure 3. (color online). Spectral functions as obtained for 
model I (Ul in the one-magnon approximation: (a) Ad{kx,ijj) 
and (b) Ac{kx,uj) ((5|. The dashed line indicates the static 
hole energy eo- Parameters; a = |, t/J = 0.25 and r] — O.lt. 



result from band folding associated with the two-site unit 
cell, independent of the ID nature of hole motion. These 
states give two QPs at each momentum in yl^*^-* (k, w), 
except a,t kx = TT where er , = eo- In contrast, the lower 



energy QP, e^^, is hidden in Ac"'' (k, w), while the higher 
energy one, £^_, has twice larger spectral weight, due to 
interference between the two sublattice contributions. 

The difference between Ad{k,uj) and Ac{k,uj) (O for 
model I which follows from states parity is well visible 
in the weak coupling regime of J ^ see Fig. [3] [s^ . 
At low energy, w < eo, one finds QPs in Adik^^i^), but 
not in Acikxjio). The width of the QP band is some- 
what reduced from the expected 2t, but even stronger 
renormalization (its width is lower than t) is found for 
w > eo. New incoherent features observed in both spec- 
tral functions ([5]) above uj = 2t (Fig. [3|) are generated by 
a magnon excitation due to the interchain hopping t±, 
see Fig. 2(c). These states disperse on the scale of ^ 2t 
(with energy difference between the respective maxima 
at kx = and tt being ~ 1.4t) and we recognize here the 



propagation of a holon |37l . I38| . 

Iir the intermediate coupling regime t > J, the spec- 
tral weight moves to higher energies (Fig. [J]). First, one 
recognizes distinct QPs in Ad{k,uj) [hidden in Ac(fc,aj)], 
with further decreased total width of this subband and 
the spectral weight decreasing from fc^; = to tt. This 
follows from hole scattering on magnon excitations (Fig. 
2) that become gradually more important when t/J in- 
creases. Second, the spectra Ac{k,uj) at t/J = 2.0 [Fig. 
mjb)] have three notable features: (i) the suppressed QP 
peaks at the onset of the spectrum, (ii) the upper holon 
branch, with spectral weight moving to higher energies 



Figure 4. (color online). Spectral functions as obtained for 
model I JT]) in the one-magnon approximation: (a) Ad{kx,oj), 
and (b) Ac{kx,Ljj) ([5]). The inset shows fca;-dependence of the 
frequencies uji corresponding to the four distinct features vis- 
ible in Aii{kx,ijj); the lji and LJ2 ones are hidden in Ac{kx,Ld)- 
Parameters: a — ^, t/ J — 2.0 and rj = O.lt. 



when kx decreases from tt to 0, and (iii) large spectral 
weight with weak momentum dependence at the upper 
edge of the spectrum. These higher energy features are 
expected to be strongly renormalized when the constraint 
to one-magnon excitations is lifted. In contrast, the QPs 
disappear iir the strong coupling regime t ^ J as well, 
because the destructive interference is due to parity. 

Indeed, the unbiased ED [s^ performed on a periodic 
cluster of = 24 sites [Fig. 1(a)] for model I with the 
same a — ^ yields very incoherent spectral weight dis- 
tribution in Ac(k, gj), see Fig. [5lja), in contrast to the 



2D t'J model [28|, [40[. As most important feature, the 
spectral weight is moved to high energy for the P and F 
poiirts. At the M point one observes a broad feature at 
uj/t ~ 0.3, accompanied by a shoulder at uj ~ 2.5t, and 
a small QP peak at w ~ — 2.5i — such a peak is also ob- 
served at K but absent for momenta inside the Brillouin 
zone (BZ). The spectral weight is transferred from high 
to lower energy when one moves from the F point to the 
edge of the BZ (M and K point). These qualitative fea- 
tures are generic and hold in a broad range of parameters 
(see the Appendix). 

The spectral weight is distributed quite differently in 
Ad{k,uj) [Fig. mb)]. Here QP peaks appear for all mo- 
menta at energies uj < —2t, with much more weight than 
those observed in Ac(fc,w), and there is strong spectral 
weight transfer to low energies at P and F points. The 
QP weight is maximal for the F point and small for K and 
Af, as expected from the variational approach [Fig. IH^a)]. 
We suggest that the lowest energy QP at the F point de- 
termines the chemical potential ^ and is responsible for 
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Figure 5. (color online). Spectral functions as obtained by 
ED for N = 24 sites at selected high-symmetry points F, K, 
M, and P in the BZ at t/J = 5 and r) = O.lt: (a) Ac{k, uj) and 
(b) Ad{k,oj), both for model I (1) with a = |; (c) Ac{k,oj) 
for model II (3), with a = 0.4, j2 = 0.2 and js = 0.5. ARPES 
spectra for Na2lr03 [nj with background subtracted [43 | are 
shown in (d) for the selected k points [i^ . 



the experimentally seen pseudogap of A ^ 0.35eV^ [11]. 

It is quite remarkable that the ED results obtained 
with model II [Fig. [He)] for parameters suggested [l^ 
for Na2lr03, are qualitatively similar to those found in 
model I [Fig. [5t^a)[, but have somewhat richer structure. 
At the M point one finds a two-peak structure with a 
weak QP at the low energy side, which develops to a 
shoulder at the K point. The spectra found at P and 
F points look similar to those of model I, again with 
most of the weight concentrated at high energy. The 
broad incoherent spectral weight and its shift to higher 
energy between the K and F point are recognized here 
as universal features for the parameters favoring zigzag 
phase. Indeed, moderate changes of parameters (a, and 
for model II j2 and js) modifying the degree of frustration 
while still supporting the zigzag order, do not lead to 
emerging QPs (see the Appendix). 

We suggest that the ED [39| simulates here the main 
features of the ARPES spectra measured at T = 130 
K \Tl\ — although the Neel temperature Tjv ~ 15 K is 
much smaller due to frustration, the incoherent part of 



the spectra should be only weakly affected by thermal 
fluctuations for T < J. The spectra are incoherent and 
no QPs are seen at the low edge uj ^ fi, see Fig. E^d). 
With the total width of the spectrum ~ 6t we estimate 
that t ~ 0.3 eV. Large spectral weight at the F and P 
points is seen mostly at high energy, as in the ED in 
Ac{k,uj) but not in Ad{k,uj) [Figs. [5la) and[5tb)[. At K 
and M one recognizes three characteristic features with 
appreciable spectral weight in the ARPES data [lH at 
(jj — II ~ 0.6 '411, ~ 1.1, and ~ 1.5 eV which have some 
correspondence in the ED spectra. The basic features 
of the ARPES spectra, (i) the strong incoherence and 
(ii) the shift of the spectral weight to high energies at the 
P and F points, are dictated by the zigzag correlations 
and are similar in the two models. Yet, we notice differ- 
ences in the fine-structure which reflect the dramatically 
different underlying Hamiltonians, and therefore suggest 
that model I is closer to the experimental data ll|. 

Summarizing, we have shown that essential features 
seen in recent ARPES experiments for Na2lr03 ll| may 
be described by a universal phenomenology which does 
not depend on details of modeling. The framework used 
is a t-J-like model with nearest neighbor Kitaev and 
Heisenberg exchange, having conceptually some similar- 
ity to high Tc cuprates. The physical spectral function 
Ac(fc,cij) in the strong coupling regime explains qualita- 
tively the incoherent nature of the ARPES spectra, with 
large spectral weight at high energy and the pseudogap 
determined by hidden quasiparticles at low energy. This 
also confirms that the low-energy electronic structure of 
Na2lr03 can be described by the motion of composite 
J = i fermions due to strong spin-orbit coupling. 
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APPENDIX: PARAMETER DEPENDENCE 

This Appendix presents additional numerical evidence 
which demonstrates that the results of exact diagonal- 
ization (ED) presented in Fig. 5 for the two models 
considered in the Letter are robust. In each case the 
spectra only weakly depend on the model parameters in 
the regime of the zigzag magnetic phase. 

For convenience, we repeat first the definitions of the 
models considered in the ED calculations. We consider 
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Figure 6. Spectral functions Ac{k,ijj) as obtained by ED for 
a periodic cluster of A*' = 24 sites at selected high-symmetry 
points r, K, M, and P in the BZ at t/J = 5 and 77 = O.lt for 
model I Eq. (1), with: (a) q = |, and (b) q = |. 



Figure 7. Spectral functions Ac{k,ijj) as obtained by ED for 
TV = 24 sites at selected high-symmetry points F, K, M, and 
P in the BZ at i/J = 5 and i] = O.lt for model II Eq. (3), 
with j2 = js = 0.8, and: (a) a = 0.6, (b) a = 0.4. 



the following i- J-like model /, 

ntj = nt + Hj = tY, els- 

{«i>o" 

+ J^^5^5..5,-^JkEs7'5;. (6) 

on the honeycomb lattice. The operators in the hop- 
ping term Jit, with i > 0, stand for creation of composite 
fermions with pseudospin flavor tr in the restricted space 
without double occupancies; the constraint follows from 
large on-site Coulomb repulsion U. The exchange term 
Hj describes in this case the Heisenberg-Kitaev model, 
with ferromagnetic (FM) Heisenberg {Jh < 0) and anti- 
ferromagnetic (AF) Kitaev {Jk > 0) exchange, given by 
the interaction J and a £ [0, 1], see Eqs. ([2]). The pa- 
rameter a interpolates between the Heisenberg (a = 0) 
and Kitaev (a — I) exchange couplings acting between 
pairs of nearest neighbor (NN) spins 5* = i; zigzag AF 
order is found in a broad range of a . 

The second model investigated in the Letter, called 
model II, has inverted signs of the NN spin couplings in 
the —'Hj term, and contains in addition exchange inter- 
actions between next NN (NNN) and third NN (3NN) 
spins, J„ = J(l — a)jn > with n = 2, 3 (2^: 

n'tj = Ht-Hj + [j2 S,-S, + J3 '^,•5,}. (7) 

{ijjGNNN {ii}63NN 

Further neighbor Heisenberg interactions {J2,J3} are 
necessary to stabilize zigzag AF order which is unsta- 
ble in case of only NN AF Heisenberg and FM Kitaev 
interactions (at J2 = J3 = 0) (23 |. 

The ED calculations were performed on a periodic clus- 
ter of = 24 sites for model I and model II for the pa- 
rameters favoring the expected zigzag magnetic order in 



the ground state. Note that spin correlations are system- 
atically reduced, compared to those in a zigzag-ordered 
phase, by the fact that no symmetry is spontaneously 
broken in the ground state of a finite cluster. First, we 
vary the parameter a and change the balance between 
the FM Heisenberg and AF Kitaev interactions in Eqs. 

from a — ^ to a — ^. As the energy increases in this 
range with increasing a [24], frustration of spin interac- 
tions in Eq. ([1]) increases as well. Surprisingly, the spec- 
tra presented in Fig. 6 and in Fig. 5(a) for the selected k 
points in the Brillouin zone (BZ) are only rather weakly 
influenced by that. The overall width of somewhat less 
than 6t does not change showing that it is determined 
by the hopping term Hf The spectral weight is mostly 
incoherent — it concentrates in the central intermediate 
energy part for the M and K points at the edge of the 
BZ, while it is moved to high energy for the P and F 
points. A weak low energy quasiparticle (QP) is visible 
at the M and K points with the spectral weight growing 
with a. But even at large a = | these QPs are rather 
weak and incoherent spectral weight dominates. The two 
cases displayed in Fig. 6 are nearly identical to each other 
and to the representative spectrum for a = | presented 
in Fig. 5(a). 

Second, we consider model II Eq. ([7]) and remark that 
the ED spectra obtained with it are qualitatively similar 
to those of model I at the same points in the BZ, see Fig. 
7. Again, the spectral weight is mostly incoherent and 
moved to high energy at the F and P points. Here we 
consider two interaction parameter sets {a, j2, Js} which 
have less frustration than the case of a = 0.4, j2 = 0.2 
and j3 = 0, 5 shown in Fig. 5(c), and the zigzag magnetic 
order is even more favored. While the spectral functions 
at the F and P points are nearly insensitive to these pa- 
rameter changes, a weak QP forms and is better resolved 
at the M point for both parameter sets of Fig. 7 than 
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at {a,j2,j3} = {0.4,0.2,0.5}. Again, the spectra found 
at the points belonging to the edge of the BZ {K and 
M) have richer structure than those in model I (Fig. 6). 
At the M point a two-peak structure is found, with a 
weak QP forming a shoulder of the lower peak at the low 
energy side. The broad incoherent spectral weight and 
its shift to higher energy between the edge and the cen- 
ter of the BZ are recognized as universal features for the 
parameters favoring zigzag magnetic order. 

Summarizing, the numerical data show that for both 
models the suppression of low energy QPs occurs in a 
broad range of parameters which correspond to the zigzag 
magnetic order in the ground state, and is therefore not 
a result of parameter tuning. This makes the observa- 
tions made in the Letter robust features of the presented 
models I and II. We have found: (i) the incoherent dis- 
tribution of the spectral weight, (ii) large spectral weight 
at high energy, and (iii) the absence of low energy QPs 
from the ARPES spectra because of destructive interfer- 
ence, even though QP eigenstates exist in the low-energy 
spectrum. This could explain pseudogap-type behavior 
in these materials. 

These results suggest that the qualitative agreement 
with the recently reported experimental angle-resolved 
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photoemission spectroscopy (ARPES) spectra 
generic, and thus bring new insights into spectral prop- 
erties of such a strong spin-orbit-coupled Mott insulator. 
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